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EFFECT OF STREAMWISE VORTICES ON 
TOLLMI EN-SCHL ICHTING WAVES 
Ali H. Nay f eh 

Department of Engineerinn r.ience and Mechanics 
Virginia Polytechnic Instil \e and State University 
Blacksburg, Virginia 24061 

The method of multiple scales is used to determine a first-order 
uniform expansion for the effect of counter-rotating steady streamwlse 
vortices in growing boundary layers on Tollmien-Schlichting waves. The 
results show that such vortices have a strong tendency to amplify 
three-dimensional Tollmien-Schlichting waves having a spanwise wave- 
length that is twice the wavelength of the vortices. An analytical 
expression is derived for the growth rates of these waves. These growth 
rates increase linearly with increasing amplitudes of the vortices. 

I. INTRODUCTION 

We consider the effect of counter-rotating steady streamwise 
vortices on the instability of growing boundary layers. We describe 
a parametric instability mechanism by which such vorticies amplify 
selected three-dimensional Tollmien-Schlichting waves. To first order, 
the selected waves have a spanwise wavelength that is twice that of the 
vortices. 

Weak and moderately strong steady streamwise vortices arise abun- 
dantly in boundary layers from many causes. In a series of wind-tunnel 
tests over flat plates with zero-pressure gradient, Klebanoff & Tidstrom 
(1959) observed steady quasi-^eriodic variations in the spanwise direc- 
tion (streamwise vortices) evidently evoked by freestream conditions. 
Similar vortices were observed in a National Physical Laboratory tunnel 
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specifically designed for the study of two-dimensional boundary layers. 
Bradshaw (1965) found that these variations may appear downstream of 
slightly nonuniform settling-chamber damping screens, depending on their 
solidity. Using the method of matched asymptotic expansions, Crow (1966) 
inferred the effect of a small, periodic incident transverse flow on the 
mean boundary layer over a flat plate. 

Gortler (1941) found that a boundary layer over a concave surface 
is strongly unstable. The instability is manifested by the presence of 
counter-rotating vortices having their axis in the streamwise direction. 

Using the tellurium method, Wortman (1969) gave detailed flow visualization 
of these vortices on slightly curved walls. Bippes (1978) conducted 
experiments on walls with radii of curvature of 0.5 and lm so that the 
generated Gb'rtler vortices are fairly strong. He made the flow visible 
by using the hydrogen-bubble technique and photographed it with a photo- 
grogrammatic sterocamera. He analyzed photogrammetrically the photographs 
and obtained fairly accurate quantitative information of the flow field. 

Unlike the case of preexisting streamwise vortices, Gortler vortices generated 
by a concave surface are amplified with streamwise distance. Their 
amplification is exponential (Smith, 1955) when they are weak and it 
appears to be linear when they are strong (Bippes, 1978). Floryan and 
Saric (1979) gave a comprehensive review of the different analyses of 
these vortices and, using the modtl of Smith (1955), presented fairly 
accurate numerical results describing these vortices. 

In their experiments on the nature of boundary- layer stability, 

Klebanoff, Tidstrom & Sargent (1962) established that streamwise vortices 
are associated with nonlinear three-dimensional wave motions. Benney & 

Lin (1966) modeled the generation of these vortices by the nonlinear 


2 


interaction of a two-dimensional wave with a three-dimensional wave 
superimposed on the laminar profile. Using this model, Antar & Collins 
(1975) calculated these vortices in a boundary layer on a flat plate 
for different amplitudes of the waves. 

A number of experimental studies investigated the influence of 
steady streamwise vortices on the transition from laminar to turbulent 
flow. Aihara (1962), Tani & Sakagami (1964), and Tani & Aihara (1969) 
studied the influence of steady streamwise vortices on two-dimensional 
Tollmien-Schl ichting waves generated by a vibrating ribbon. They con- 
sidered the case in which the vortices were generated naturally on a 
concave surface (Gb'rtler vortices) as well as the case in which the 
vortices were generated artifically by a row of wings on a cancave 
surface. They measured the distributions of the mean velocity and wave 
intensity across the boundary layer for three spanwise positions at a 
number of streamwise stations. They concluded that the GSrtler vortices 
indirectly affect the transition by inducing a spanwise variation in 
boundary-layer thickness, at least when the radii of curvature are not 
extremely small. However, they did not present any measurements of 
the growth rates of the Tollmien-Schl ichting waves. Wortman (1969) 
investigated the development of natural transition downstream of Gortler 
vortices. Using the tellurium method, he determined the direction and 
relative magnitude of the unsteady velocities from the streaklines by 
confining his observations to the vicinity of the starting point of the 
streaklines. He observed a steady second-order instability that destroys 
the symmetry of the Gortler vortices. He suggested that this instability 
is caused by secondary vortices having spanwise wavelengths that are 
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twice those of the Gb'rtler vortices. Then he observed a third-order 
instability, consisting of regular three-dimensional oscillations. 

The above shows that there are many theoretical and experimental 
studies relating to the generation of streamwlse vortices and a number 
of experimental studies relating to their effect on transition, but to 
the author's knowledge, no theory yet exists on how these vortices 
affect the development of Tollmien-Schlichting waves. The purpose of 
the present paper is to present a parameteric instability mechanism by 
which the streamwise vortices increase the growth of selected Tollmien- 
Schlichting waves in growing boundary layers. To first order the 
selected waves have a spanwise wavelength that is twic* *hat of the 
vortices, while to second order the selected waves have spanwise wave- 
lengths that are equal and twice that of the vortices. To minimize the 
algebra, we consider flows with growing steady counter-rotating vortices 
over flat plates with or without pressure gradients. 
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II. Problem Formulation 


We consider the stability of a basic flow that consists of the 
superposition of the Blasius or Falkner-Skar. flow and a flow corres- 
ponding to growing steady quasi-periodic counter-rotating streamwise 
vortices. Thus we consider the stability of the flow described by 

U * UoUi.y) + EyUi (xi ,zi ,y)cos2Bz + .... 

(1) 

V = eV 0 (xi,y) + e v Vi(x 1? z ,y)cos28z + 

(2) 

W = e y Wi (xi ,Zi ,y)sin28z + ..., 

(3) 

P = Po(xi) + e v Pi(xi,Zi,y)cos2Bz + .... 

(4) 

where the subscript o refers to the Blasius of 

Falkner-Skan flow, the 


subscript 1 refers to the flow corresponding the the streamwise vortices, 
Xi = ex, Zi = ez, e is a small dimensonless quantity that accounts for 
the quasi -periodicity of the vortices as well as for the streamwise growth 
of the vortices and the boundary layer, 8 is a real dimensionless spanwise 
wavenumber, and e y is a small dimensionless quantity that indicates the 
strength of the streamwise vortices. In Eqs. (1) - (4), velocities and 
lengths are made dimensionless by using a reference velocity U r and a 
reference length 6 r - We superpose the small unsteady perturbation 
quantities e T u(x,y,z,t) , e T v(x,y,z,t) , £jw(x,y,z,t) , and e T p(x,y,z,t) 
on those given in Eqs. (1) - (4) so that the total flow quantities 
become U + e-j-u, V + CyV, W + CjW, and P + e^p. Here, Cj is a small 
dimensionless quantity that is the order of the amplitude of the Tollmien- 
Schlichting waves. In this paper, e-j- is assumed to be much smaller than 
e v and e so that terms the order of can be neglected compared with CjC y 
and Eye. Substituting these total flow quantities into the dimension- 
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less Navier-stokes equations, subtracting the basic-flow quantities, and 
keeping linear terms in the perturbation quantities, we obtain 

Uu 1 * . 

3x 3y a 2 " u » ( 5 ) 

^ +U n + “l^ + ^**|tH|i + w|U..|£ t l v2ui(6) 


If + u |i + u 2L ♦ V |* ♦ v 3V 

3t 8x 3x 3y Sy ■ " 51 • " 3Y - 57 ■ jf 


& + u 

at u 


3W 

ax 


+ u f? + '|*'f*«|?*w|f.-|| t l v!w , (8) 

where t is made dimensionless by using V U r , R = U(V v is the Reynolds 
number, and v is the kinematic viscosity of the fluid. 

Substituting Eqs. (l)-(4) into Eqs. (6)-(8). we obtain 

I? * “• & * * f 1 * |f - p V‘u = -c„ jiu, |f ♦ V, 0 ♦ ||l. ,] C os26z 
* CW ‘ I? ' 26U -'* i "2Bt] - s[||k u ♦ V, |H] ♦ 0( ee , (9) 


It * Ut & * f? - i »** * - c y ^U, |f + V, |f * !£. v]cos23z 

* t H i fj - 28V,w]sin2Bz( - e[Vo |f ♦ |i 


J 


a + U, |£ + |£ . 1 v=w • - c 
3t ax a z r w £ v 


V^* 3x 3y T 37 

v] + 0(e 2 ,ee ) 


r 


+ # v + W > fr]sin2Bz(- e V 0 & ♦ 0(eej 


00 ) 


|[ u i |~ + V 2 0 + 2 BW!w]cos2Bz 


3w 


L 


( 11 ) 


ay v- 

Equations (5) and (9)-(ll) need to be supplemented by initial and bound- 
ary conditions. The initial conditions are specified later, while the 
boundary conditions for an impermeable flat surface are 

u = v = w = 0 at y = 0, (12) 


u. V, w 


at y 


(13) 
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III. Solution 


We use the method of multiple scales (e.g., Nayfeh, 1973) to 
determine a first-order uniform expansion for Eqs. (5) and (9)-(13). To 
accomplish this, we let e y = 0(e) and write e y = where x = 0(1). If 
e << e v , the effect of the growth of the boundary layer is small com- 
pared with the effect of the vortices. If c y << e, the effect of the 
vortices is small compared with that due to the growth of the boundary 
layer, and the solution accounts for the nonparallel effects only. Thus, 
the above ordering yields an expansion that accounts for the effects of 
the streamwise vortices and the growth of the boundary layer, and it 
includes the cases e << e y and c y « e as special cases. 

We seek a uniform expansion for Eqs. (5) and (9)-(13) in the form 


u= { e n u n (x 0 ,Xi,.v,Zo,Zi.to,ti) + 0(e 2 ), 

n=0 n 


v = f e n v (x 0 ,Xi ,y,z 0 ,Zi ,t 0 ,ti ) + 0 (e 2 ), 
n=0 r 

i n 

W = l c w (xo.Xx.y.Zo.Zl.to.ti) + 0(e 2 ) , 

_ A ' * 


P = 


i « 

y e P n (x 0 ,Xi,y,z 0 ,Zi,to,ti) 
n=0 n 


+ 0(e 2 ), 


where 


x n - c x. 


z n = 


,n, 

L Z | 


t n = 


n + 
e t 


(14) 

(15) 

(16) 

(17) 

(18) 


Substituting Eqs. (14)-(18) into Eqs. (5) and (9)-(13) and equating 
coefficients of like powers e, we obtain 
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09 ) 


Order c f 


/ 2 (u 0 ,v 0 ,Po) « + Uo 1^- + Vo |^L + Ifj- - J V5U 0 - 0, 

/3(vo,Po) = llj- + u o fj£ + V 2 oV 0 « 0, 

/4( w 0,Po) « !§■ + Uo IS' + - i V 2 0*0 * 0, 


Uo = V 0 = Wo = 0 


Uo • Vo ,w 0 "* 0 


3x 0 ’ 3Zo 
at y = 0, 

as y ■+■ », 


(20) 

( 21 ) 

(22) 

(23) 

(24) 


Order e 


/,(u,, v ,.w,) * - if -ijf . 


• xjiu, if ♦ Vi §*■ * i 1 - V,]C0S26Z, * [W, 
- ZSU^cJs'n^SZoJ - 


9Uq 

3z 0 


3Uo 
3x i 


Uo - Vo 


(25) 


(26) 


3v 


^ 3(Vl,Pl) = • It' Uo Ixf + S c I^3“x 7 + liTaTT 3 ' x f [Ul 9Xo 

*“■ -) 

+ V , |^2- + v 0 ]cos2Bz 0 + [W, - 2BViW 0 ]sin2Bz 0 

- v °!r-f Lvo * 

r 

[ y i $5j' + Vl + 26 WiW 0 ]cos2Bz 0 + [|^ v 0 


(27) 


+ W, f"]sin 26 zj - 


J 


3 W[ 

3 y 


(28) 
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Uj * Vj * wt * 0 at y ■ 0, 
ui.vj.wi 0 as y ■*> », 


(29) 

(30) 


where 


v2 ° “ m * w + Hr • 


The initial conditions are taken such that the solution of the 
zeroth-order problem, Eqs. (19)-(24), consists of two wave packets 
centered around the frequency w, the streamwlse wavenumber a, and the 
spanwlse wavenumbers Si and - Si; that 1$, 

Uo * Aj (xj ,z i ,t i )c 1 1 (y ,Xj )exp( 19 1 ) + Aj(xj »Z} * tj )^ 2 1 (y*Xj )exp(i 02 ) * 

(31) 

vo r A i c i 2 (y» x i )exp(i6i ) ♦ A 2 c 2 2 (y,xi)exp(ie 2 ), (32) 


W 0 * A,r,i 3 (y.xi)exp(16i) + A 2 c 2) (y,Xi)exp(i6 2 ), 


Po s AiCm(y,Xi)exp(i6i) + A 2 c 21 ,(y,xi)exp(ie 2 ), 


where 

6 n -y 


G„ ■ / ot(x j )dx - S n z 0 - wto, S 2 * - Si , 


(33) 


(34) 


(35) 


and the functions Ai and A 2 are undetermined at this level of approxima- 
tion; they are determined by imposing the solvability conditions at the 
next level of approximation. Substituting Eqs. (31 )- ( 35) into Eqs. 

(19)- (24) yields the following eigenvalue problems: 
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( 36 ) 


1o ' n . + <*„, • < 8 n c n , * 0 • 

1(U,o • ♦ c ni ou. ♦ ioC ni - J (O 1 - o‘ - ei)Cn, ■ 0. (3?) 

1(U,0 - u); n! ♦ Dc n< - J (o 1 - o' - ei)e ni * o. (38) 

HU go - u )c n> - 18„; n< - if 0> ! - o ! - 6()C ni ■ 0. (39) 

«!., ■ ■ «n, ■ 0 ** » ■ °- <40) 

VV'n. * 0 (,,) 

where Dc - 3r./3y. For a given w. Si, and k, one can solve Eqs. (36)- 
(41) numerically to determine the complex eigenvalue a and the eigen- 
functions 

nm 

Substituting Eqs. (31)-(35) Into Eqs. (25)-(28) yields 

(“>.*>•“») ■ - L *3Sr«n, 

n* I 

- Jl (42) 

3 A 3A 

c ^2(ui,v l ,p l ) * - i ; ni ♦ (U«; ni + ’ ir ; n^ 3x7 

n* I 

* 2 Vn, 57"5 ex p()8 n ) - j t C(U, - Hr> * R 377 ] 

* 5^*- * it Si, * v * ir 11 V ,p(,9 n ) ' 1 xC'U.o;.i 

♦ v, fixx * - W,6,c„ ♦ zmi.o.w.ntfi*. 

♦ 21(3 - 6j)Zo] - \ xC^UjOCj! ♦ Vi C 22 

- WjBiC ? j - 2iBU,r, 2J ]A 2 exp[1e 1 - 21(6 - BiUo] ♦ NST, (43) 
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X 


e *, (, ll p,) . . C n; ♦ (u, - ^H ni 

«xp(ie ) - f c(u. • ir) sr^ - 5 c 5 s - 
" n »1 R 3x » * n? dxT 

+ 3y ^ V ° r, n 7 ^ exp ( 1 V ” 7 x[*«UiCia + |p“( V i c a 2 ) - gtWiCu 

+ 216ViCi j]Aiexp[1e 2 ♦ 21(8 * e,) 2 0 ] - \ X [1aU 1 C 22 

+ ay ( v iC 2 ?) - SiWiCz? - 21 BV 1 ^7 3 3 A?expCl e 1 - 21(8 - 81)20] 


+ NST, 


(44) 


3A ^ 

«r 1 “ p(,# » ) ; J, [(U ° • TT> 5^ - * 37T 

+ Vo 3^]A n exp(1e n ) - ^rxClUiad, ♦ (28 - 8i)W,eu 


+ V, 


fy^ ' 1 f - 1 Ci 2 ]Aiexp[i 0 7 + 21(8 - 81 ) 20 ] 


* l x[iU,ar. 7J + (28 - 8 i)W l ; ? , + V, 

3 W 

* 1 ^ C??]A ? exp[ 10 i - 21(8 - 61 ) 20 ], 


r ’n 3 ^ 


(45) 


where NST stands for terms that are proportional to exp[?i (6 + 8 i)z 0 ], 
which do not produce secular terms in uj, Vi, w Jt and p,. 

Since the homogeneous parts of Eqs. (42)-(45), (29), and (30) are 
the same as Eqs. (19)-(24) and since the latter have a nontrivial solu- 
tion, the inhomogeneous equations (42)-(45), (29), and (30) have a 
solution if, and only if, the inhomogeneous pares are orthogonal to 
every solution of the adjoint homogeneous problem, gjhese solvaH'..ty 
conditions depend on whether 8 ~ 61 or not. If 6 is away from e lt 
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the solvability conditions yield two uncoupled equations describing the 
effect of nonparallelism on Aj and A 2 . If B : Bi, we Introduce a 
detuning parameter o defined by 

6 = Bi + eo ( 46 ) 

where a 0(1) and express (B - 2Bi)z 0 as azi. Then, imposing the 
solvability condition that the inhomogeneities be orthogonal to every 
solution of the adjoint homogeneous problem, we obtain 

911 fir + 912 IxT + 913 IzT = hlAl + xhi 2 A 2 exp(-iazi) (47) 

921 3tf + 922 1x7 + 923 flf = hzA 2 + Xh 2 iAiexp(iazi) (48) 

where the g's and h's are given in Appendix A together with the adjoint 
problems. Differentiating Eqs. (19) - (24), respectively, with respect 
to a and B n and imposing the solvability conditions, one can show that 


_n» 

< ‘ ““ 


fru _ 

9 nj " \ 


(49) 


where and are the complex group velocities in the x and z directions 
Since the solutions of Eqs. (47) and (48) for general initial 
conditions are not available yet, we consider next the special case of 
periodic streamwise vortices and a single-frequency disturbance that is 
perfectly tuned in the spanwise wavenumber. The single-frequency 
assumption corresponds to the case of a disturbance generated by a 
vibrating ribbon. The second assumption demands that B = Bi and that 
the waves are modulated in the streamwise position only. Thus, we 
consider the case in which 3A n /3t, = 3A n /3 Zl = 0 and o * 0. Then, Eqs. 

(47) and (48) can be rewritten as 
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( 50 ) 


* ehiAi + e y hi 2 A 2 
gj 2 -- eh 2 A 2 + e v h 21 Ai 


(51) 


where 


h n = h n /g n2> hj 2 = hi 2 /gi 2 , h 2 j = h 2 i/g 22 

It follows from Eqs. (36) - (41) that Ci 3 = - c 23 and 

4i n s C 2 n for n ^ 3, while it follows from Eqs. (A. 7) - (A. 12) that 

★ ★ ★ ★ 

the adjoint solutions are related by c 14 = - c 24 and Ci n = c 2n for 
n f 4. Hence, it follows from Eqs. (A. 2), (A. 4), (A. 5), and (A. 6) that 

9 1 2 * 922 . hi = h 2 , hi 2 = h 21 . 

Thus, 

A A A A 

hi = h 2 and hi 2 = h 2 i» 

Therefore, adding Eqs. (50) and (51) yields 


(A 2 + Ai) = (eh 1 + e y hi 2 )(A 2 + Ai) 

(52) 

Subtracting Eq. (50) from Eq. (51) yields 


g^- (A 2 - Ai) = (chi - e y hi 2 )(A 2 - Ai ) 

(53) 

The solutions of Eqs. (52) and (53) are 


M As As 

A 2 + Ai = 2ciexp[ J (chi + e y hi 2 )dx] 

(54) 

r a a 

A 2 - Aj = 2c 2 exp[ J (chi - e y hi 2 )dx] 

(55) 


where Ci and c 2 are arbitrary constants that can be determined from the 
initial conditions. Solving Eqs. (54) and (55) gives 

/ A A 0 A A 

(eh x + e v hi 2 )dx] + c 2 exp[ J (ehj - e y hi 2 )dx] 

Ai = Cjexp[ J* (chi + e y hi 2 )dx] - c 2 exp[ J (chi - c y h J 2 )dx] 

(57) 
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Substituting for Ai and A 2 in Eq. (31), using Eq. (35), substituting 
the results into Eq. (14), and recalling that Bi * 8, we obtain 


u =£nexp[ /<- + ehi)dx - iez - iwt] x 
[ciexp(e v J*h X2 dx) - c 2 exp(-e y J*hi 2 dx)] 


+ C21 


exp[ J*l 


(ia + ehjdx + i$z - iait] x 


[ciexp(e v J*hi 2 dx) + c 2 exp(-e y J*hi2dx)] + ... 


(58) 


Equations (A. 2) and (A. 5) show that, in general, h X2 is a complex 
number. Hence, one of the terms multipling ci and c 2 decays while the 
other grows exponentially with distance. Thus, the growth rate 0 based 
on u for either the wave with the positive or negative 6 is 


a = - ^ + e[Real(hi) + f^] + e y | Real(hi 2 )| (59) 


because Cu = 4 2 i . 

Equation (59) shows that the growth rate is the sum of three 
quantities: - , the quasi -paral lei growth rate; the term proportional 

to e, the effect of nonparallelism; and the term proportional to e y , the 
effect of the streamwise vortices. Thus, in a given physical situation, 
the relative influence of the vortices and nonparallelism depends on 
the relative magnitudes of e and e y . For maximum amplified waves, 
e = 0(10 ), while for flows over concave surfaces, e y can be (0.10), 

depending on the radius of curvature. In such situations, the effects of 
the vortices dominate the effects of nonparallelism, and the presence 
of the vortices is a very powerful instability mechanism. Hence, the 
presence of this mechanism may not be difficult to check experimentally. 
This requires an experiment in which the amplitudes and spanwise varia- 
tions of Gortler vortices and Tol lmien-Schl ichting waves on a concave 
wall are measured. The spanwise variations of the Gb'rtler vortices and 
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the generated Tollmien-Schlichting waves can be checked to see whether 
the wavelengths of the latter are twice those of the vortices. The 
modification of the mean flow due to the presence of the vortices can be 
determined and the interaction coefficient h 12 can be calculated from 
Eqs. (A. 2) and (A. 5) by quadrature. Then, the variation of u with x and 
y can be calculated from Eq. (58) and compared with the experimental 
results. 

It should be noted that the present analysis is valid only when 
the amplitude of the Tollmien-Schlichting waves Cy is small compared 
with the amplitude of the vortices e y . As the Tollmien-Schlichting 
waves grow, one needs to account for their influence on the vortices. 

In fact, they will generate streamwise vortices having an amplitude 
0(e|) (Klebanoff, Tidstrom & Sargent, 1962; Benney & Lin, 1960), which 
may strengthen or weeken the primary vortices, depending on their 
phasings. This effect has not been taken into account in this paper. 
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APPENDIX A 


'rh ~ ^^1 n! C n 2 + ? n 2 ^n 3 + ? n 3 ^n 
«/o 

n, ■ J tvs. * + (“» - rut,,';, + vs, ♦ v;, )]<iy 

> 210 
n 3 "J ^n 3 c n! + c n 4 c n 4 + ~R ^ 


(A. 2) 

ni ? n 2 + ^n 2 ^n 3 + ? n 3 5 n 4 ^ dy ^ A * 3 ^ 


3 t r 


W A 

- h =f[!!nir* + [(u. - ila.) !5ai.- i, ^ + 

n J J8X! Sj u 0 R 9 xj R c n 3 dxj 3 x 3 3 x x ? n x 
+ V ° ay^K + Wo ' TT^ 3x7^ ' f c n 2 Bx7 + |y < V « V^ 3 


2 ia x 9 c n 3 i „ da 9c 


- 2 h 


* 2 h 


* [(u " - ir> 5^ - i % ir + Vt ^S,] dy « A) 

‘ °1 ^’“ Ul!!1 * V ‘ ly^* fy'-t” ' B ‘ w ‘ 5 *‘ ■ 21 0 u iC 2 jDcTj 

+ CiaUx ; 22 + Ij- (VjC 22 ) - S 1 W l522 - 2 i 6 V 1 ; 23 ]^3 

+ [ i aU i c 2 3 + (23 - @i )WiC 23 + v i + i S2 2 ]chjdy (A. 5) 
00 r 

. = rkiotUi?,! + Vi |iu- + It 2 - Ci 2 - eiMitn + 21 eUiCi 3 ]cJ 2 
+ [iaUiCi 2 + It (Vi^! 2 ) - $iWiCi 2 + 2 l 0 V 1 Gi 3 ]ef 3 + [iaUi?i 3 


+ (20 - 0i)W lC i, + Vi ^ - i c 12 kl v d y 


(A. 6) 
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where the c* are solutions of the following adjoint problems: 
nm 

iaC* t + i(Uoo - w)C n2 - ^ (D 2 - a 2 - 0jk n2 = 0 (A. 7) 

- D C + DUo C + i(Uoa ' “ k n 3 ' l ( ° 2 ‘ a2 ‘ B J k n3 = 0 


(A. 8) 


" ie n c ni + i(Uoa " “^n 4 " l ^ ‘ a2 " B n )c n 4 = 0 (A ' 9) 


1a C " D? n 3 * i6 "^ 


rrn 4 


c n 2 = ^n 3 


= C, 


n 4 


at y = 0 


’nm 


+ 0 as 


(A. 10) 
(A.ll) 

(A. 12) 
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